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O ' We present measurements of the [On] emission line width for a com- 

^ \ plete sample of 24 blue field galaxies (21.25 < B < 22, B - R < 1.2) 

a ■ 

at (z) ~ 0.25, obtained with the AUTOFIB fibre spectrograph on the 
Anglo- Australian Telescope (A AT). Most emission lines are spectrally 
$H ' resolved, yet all have dispersions a v < 100 km s _1 . Five of the 24 sam- 



ple members have [On] doublet line flux ratios which imply gas den- 
sities in excess of 100 cm" 3 . The line emission in these galaxies may 
be dominated by an active nucleus and the galaxies have been elimi- 
nated from the subsequent analysis. The remaining 19 linewidths are 
too large by a factor of two (7a significance) to be attributed to tur- 
bulent motions within an individual star forming region, and therefore 
most likely reflect the orbital motion of ionized gas in the galaxy. We 
use Fabry-Perot observations of nearby galaxies to construct simulated 
datasets that mimic our observational setup at z ~ 0.25; these allow us 
to compute the expected distribution of (observable) linewidths a v for 
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a galaxy of a given "true" (optical) rotation speed v c . These simulations 
include the effects of random viewing angles, clumpy line emission, fi- 
nite fibre aperture, and internal dust extinction on the emission line 
profile. We assume a lincwidth- luminosity colour relation: 



log 



v c (M B , B-R) = Iog[« c (-19, 1)] - r, (M B + 19) + C \{B-R)-1] 



and determine the range of parameters consistent with our data. We 
find a mean rotation speed of v c (— 19, 1) = 66 ± 8 km s _1 (68% confi- 
dence limits) for the distant galaxies with Mb = — 19 and B — R = 1, 
with a magnitude dependence for v c of rj — 0.07 ± 0.08, and a colour 
dependence of £ = 0.28 ± 0.25. Through comparison with several lo- 
cal samples we show that this value of v c {— 19, 1) is significantly lower 
than the optical rotation speed of present-day galaxies with the same 
absolute magnitude and rest frame colour (w 105 km s _1 ). The most 
straightforward interpretation is that the distant blue, sub-L* galaxies 
are about 1.5 mag brighter (and > 0.8 mag brighter at 99% confidence) 
than local galaxies of the same linewidth and colour. 
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1 INTRODUCTION 

In the last decade, large ground-based telescopes, HST and sensitive detectors have 
made it possible to study galaxy properties over cosmological distances in order to 
obtain direct constraints on galaxy evolution. One of the most surprising results of 
these studies has been the high number counts of "blue" field galaxies at magnitudes 
fainter than B ~ 20 (Kron 1980; Tyson 1988). Even though these counts suggest 
significant evolutionary effects, the shape of the measured redshift distribution of 
blue-selected galaxies (B < 23) is consistent with no-evolution models (Broadhurst 
et al. 1988; Colless et al. 1990, hereafter referred to as CETH; Lilly et al. 1991; Colless 
et al. 1993b). Further, deep imaging surveys in the red or near infrared bandpasses 
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yield galaxy counts which do not indicate strong evolutionary effects (Cowie et al. 
f993; Gardner et al. 1993). 

These puzzling observational results have prompted a large number of models to 
explain the observations. These models fall into two broad classes in their approach. 
First, there are ad hoc models, which assume some evolution of galaxy properties 
(number density, luminosity, etc.) and test which of these scenarios are consistent with 
the observations. The evolutionary effects in these models are not directly related to 
the physics of galaxy formation. The no-evolution hypothesis, mild and strong lumi- 
nosity and spectral evolution models (Tinsley 1972; Guiderdoni & Rocca-Volmerange 
1990; Gronwall & Koo 1995), the differential luminosity evolution model (Broadhurst 
et al. 1988), the rapidly fading/disappearing dwarf galaxy picture (Babul & Rees 
1991), and number density evolution models (cf. the merging scenario of Broadhurst 
et al. 1992) fall into the first category. A second class of evolutionary models tries to 
incorporate a substantial amount of independent information about the hierarchical 
formation of structure in the Universe (cf. Kauffmann et al. 1994; Cole et al. 1994). 
There the number density of collapsed, virialized structures is derived as a function 
of epoch from a given cosmogonic scenario (e.g., a cold dark matter cosmogony), and 
combined with prescriptions for the merging and star-formation history of galaxies. 
The luminosity and density evolution of galaxies are coupled in this latter class of 
models. 

Even though these models differ widely in many aspects, they share a common 
feature: the bulk of the evolutionary changes are attributed to intrinsically faint 
(L < 0.5 L*), blue galaxies. This feature can explain the paucity of galaxies with 
z ;> 1 in redshift surveys of 5-selected galaxies, and may explain why evolutionary 
effects appear to be stronger in observations at shorter wavelengths. The observed 
number counts and redshift distribution constrain the overall evolution of the galaxy 
luminosity function, but do not provide a unique prescription for the evolutionary 
history of individual galaxies. Therefore it is still a wide open question to which 
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extent the evolution of the luminosity function is caused by changes in the luminosity 
of individual galaxies or by changes in number density of visible galaxies. 

The key to understanding the evolution of individual galaxies lies in the identifi- 
cation of the local counterparts of distant galaxies. In this paper we try to establish 
a link between distant and local galaxies by comparing the relation between the 
spatially integrated Doppler linewidth (LW) of the ionized gas to the continuum lu- 
minosity in the two populations. This observational test effectively seeks to determine 
how the luminosity of galaxies of a given LW or mass scale — i.e., the specific lumi- 
nosity of galaxies — changes from moderate redshifts to the present day and thereby 
addresses the question of "luminosity evolution" directly. 

Various workers have shown that it is feasible to obtain line-width measure- 
ments and rotation curves for galaxies at cosmologically significant distances (Vogt et 
al. 1993, Colless 1994, Koo et al. 1995; Rix, Colless and Guhathakurta, 1995; Guzman 
et al. 1996; Vogt et al. 1996). The are three important new elements in this paper: 
(1) a statistically well defined sample, (2) detailed modeling of the observations, and 
(3) a careful comparison with properly chosen local galaxy samples. 

We have carried out an experiment to determine the global LW of the ionized 
gas — using the [On] doublet — in a complete, magnitude-limited sample of blue field 
galaxies at (z) ~ 0.25. Using a fibre spectrograph, we have performed the optical 
analog of a single beam Hi LW measurement. The target galaxies have the magni- 
tudes and colours of local, small, star-forming spiral galaxies (similar to the Large 
Magellanic Cloud), and we compare the measured LWs of these target galaxies with 
the LWs of probable local counterparts, after correcting for several biases inherent in 
such measurements. 

The paper is organized as follows: Sec. 2 contains a description of the target sam- 
ple, observations, and data reduction, and a discussion of how the LWs are measured 
from the data. In Sec. 3, we interpret the measured LW is terms of the characteristic 
rotation speed of the galaxy, and study the correlation of the galaxies' rotation speed 
with their absolute luminosity and colour. In Sec. 4, we compare the rotation speeds 
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of distant galaxies to those of local counterparts in order to quantify the evolution of 
field galaxies. The main points of this paper are summarized in Sec. 5. 

A Hubble constant of H = 70 km s^ 1 Mpc -1 and flat space-time (q = 0.5) are 
assumed throughout. 

2 SAMPLE SELECTION, OBSERVATIONS, AND DATA 
REDUCTION 

2.1 The Target Sample 

The target galaxies were chosen at random from the photometric database used by 
CETH for the LDSS-1 survey, on the basis of the following criteria: (a) apparent blue 
magnitudes in the range 21.25 < bj < 22, and (b) colours bluer than bj — r F < 1.2. 
(For the sake of convenience, we use the symbols B and R in the rest of the paper 
when referring to magnitudes in the bj and r F bands, as defined in CETH.) In 
addition to the magnitude and colour criteria used to define the target sample, the 
final sample of LWs includes only galaxies with moderate to strong [On] emission 
(equivalent widths > 10 A) within the observed redshift range 0.16 < z < 0.37, which 
corresponds to the spectral window for [On] in our setup. In Sec. 2.4, we demonstrate 
that for most galaxies that satisfy criteria (a) and (b) the requirement of detectable 
[On] emission does not impose any additional restriction. 

Our colour cut corresponds to the rest frame colours of Scd galaxies at a redshift 
of 0.25 (cf. Fig. 12 in CETH), and has been designed to yield a high fraction of galaxies 
with detectable emission lines. The photometric database from which the galaxies for 
CETH's LDSS-1 redshift survey were drawn (Jones et al. 1991) contains 764 galaxies 
in our chosen apparent magnitude range, 21.25 < B < 22, and 463 (= 61% ± 3%) of 
these are in our chosen colour range, B — R < 1.2. The colour distribution of nearly 
2000 galaxies with 21 < B < 22 in Kron's (1980) photometric sample indicates that 
<; 50% of the galaxies within this apparent magnitude range survive our colour cut. 
The median galaxy in our sample (B = 21.7 at z = 0.25) has an absolute blue 
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magnitude of M B ~ -19.0, or 0.25 L* B [L* B = -20.6 for H = 70 km s _1 Mpc" 1 
(Efstathiou et al. 1988)]. The absolute magnitudes of the galaxies listed in Table 1 
were calculated from their apparent B magnitudes, accounting for their luminosity 
distance and a K-correction for Scd galaxies (cf. Fig. 13.7 in Peebles 1993). 

2.2 Observations 

The high-resolution spectra presented in this paper were obtained with the AUT- 
OFIB multifibre spectrograph at the Anglo-Australian Telescope (AAT) on 1993 
October 12, with seeing of ^ 1 '.'7. A 1200 lines mm -1 grating was used to cover 
the spectral region from 4300 A to 5110 A at a dispersion of 0.79 A pixel -1 . The 
wavelength coverage was chosen to detect redshifted [On] (3727 A) emission line dou- 
blets in galaxies with 0.16 < z < 0.37, a redshift range centred on the peak of the 
redshift distribution of B — 21.25-22 galaxies to maximize the [On] detection rate. 
The instrumental resolution is a = 0.82 A (FWHM = 1.9 A), which corresponds to 
a = 56 km s -1 and 48 km s -1 at the blue and red ends of the spectral range, re- 
spectively. The doublet nature of the [On] emission line, with a doublet separation 
of AA = 2.75 A, is resolved at the instrumental resolution. This makes it necessary 
to tailor the data reduction process (Sec. 2.5), but has the advantage of allowing an 
unambiguous redshift determination based on only one emission line. 

The spectra were recorded on a 1024 x 1024 Tektronix CCD. The AUTOFIB 
spectrograph contains 60 usable fibres, each 2'.'1 in diameter, which can be positioned 
over a 40' diameter circular field. Engineering tests indicate that the fibres can be 
positioned with an accuracy of 0'/3 rms. The overall throughput of the system is a 
little under 1%. Of the 60 fibres, 49 were positioned on target galaxies in the FBG1 
field. The approximate coordinates of the field center are: ^1950 = 00 h 54 m 51 s ; £1950 = 
-27°51'. 

Five of the remaining fibres were positioned at random 'blank' sky locations in 
order to measure (and to subsequently subtract) the night sky spectrum, while six 
were pointed at nearby bright stars that served as astrometric references. The total 
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exposure time on the FBG1 field was 7 hr, split into individual 1 hr exposures. The 
1 hr target exposures were interspersed with short (2 min) calibration exposures of 
Cu-Ar and Fe-Ar arc lamps. The seeing was about 1 "7 (FWHM) as measured off the 
astrometric reference stars which were observed through fibre bundles. Bias frames 
and short flat field exposures of the twilight sky and of a tungsten lamp (white light 
source) were also obtained. 

2.3 Data Reduction 

The fibre spectra were reduced using standard procedures in the IRAF packages ccd- 
proc and dohydra. The individual 1 hr CCD frames were flat fielded, the wavelength 
scale in the FBG1 spectra calibrated using arc lamp spectra, and the data combined 
after rejection of cosmic ray events. A one-dimensional, optimally weighted, mean 
spectrum was extracted for each fibre from the combined CCD data, using the flat 
field image to define the response profile of each fibre and the locus of its spectrum on 
the CCD array. The spectra are dominated by night sky continuum emission outside 
the [On] lines, and the lines occupy a very small fraction (< 1%) of the full spectral 
range; the continuum flux is undetectable for most galaxies. The median of all spec- 
tra (galaxy and blank sky) was used to define an accurate template of the night sky 
spectrum and it was subtracted from each galaxy's spectrum. The sky-subtracted 
spectra of all 24 galaxies in which [On] emission is detected are shown in Figure 1. 

2.4 Detection Limits for [Oil] Emission 

2.4-1 Detection Algorithm 

We can devise an objective detection criterion for the emission lines, because galaxies 
with clearly detected [On] will not have any other emission line in the observed 
spectral range; the portion of the spectrum outside the [On] line therefore provides 
an empirical estimate of the effective noise in the data, including Poisson noise, sky 
subtraction error, residual cosmic rays, and flat fielding error. 
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We convolve the spectrum of each galaxy with a sequence of template [On] line 
profiles, each centered on the k-th pixel (with 20 < k < 1004, the portion of the 
CCD free from edge artifacts). The templates are double Gaussians of equal height 
and instrumental linewidth (cr instr = 0.82 A), broadened by various amounts to ac- 
count for the intrinsic galaxy LW. The correlation amplitude C(k), resulting from 
the convolution of these templates with the data, has a strong maximum C max at 
the redshift of any obvious [On] emission line. We calculate the rms variation, AC, 
of C(k) from the rest of the spectrum and determine the ratio of the second highest 
maximum of C(k) to AC. Since this second maximum must be due to noise, we can 
determine how high to set a threshold to avoid such spurious detections. Applying 
this statistic to all spectra with "clear" detections, i.e., with C max /AC > 10, indi- 
cates that a threshold C t h re sh/AC = 4.5 avoids any false detection. Using a narrow 
template of width, a = [of nstr + (40 km s^ 1 ) 2 ] 1 ^ 2 , we searched all 49 target galaxy 
spectra, detecting 24 galaxies with C max > C t hresh (4.5 AC). A histogram of detection 
significances is given in Figure 3 and the observed properties of these galaxies are 
listed in Table 10. 

As narrow emission lines have a higher contrast against the galaxy's continuum 
and sky background than broad lines, they are easier to detect. To check how our 
detection limit depends on LW, we have repeated the above procedure, broadening 
the template by km s _1 < o~ v < 300 km s _ . The convolution with broader templates 
resulted in no additional detections over those listed in Table 1. Further, all but one 
of the existing [On] lines could have been detected even if their LWs had been as 
large as a v < 200 km s _1 (for a given total line flux). We conclude that our detection 
efficiency is approximately uniform for a v < 200 km s~ . The fact that all of the 
24 galaxies detected have LWs a v < 100 km s" 1 is not the result of a detection bias. 

* Note that for convenience we use units of both A and km s — 1 for the dispersion cr; the ratio of these two quantities 
is A/c for the instrumental dispersion, <r; nstr , and Aq/c for the intrinsic rest frame dispersion of the [Oil] lines, a v . 
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2.4-2 Completeness 

There are two independent arguments that we have detected most, possibly all, of 
the target galaxies (21.25 < B < 22, B — R < 1.2) within the searched redshift range 
(0.16 < z < 0.37). 

(1) The strength of the emission line in most detected objects, as measured by the 
value of C max , is considerably higher than the detection threshold Cthrcsh = 4.5 AC 
The weakest detection is 5.4 AC, and the median detection amplitude is 11 AC, 
as shown in Figure 3. If most of the true [On] line fluxes for this sample of blue 
galaxies were below our detection limit and if we just saw the tip of the iceberg, 
we would have expected to find the majority of detected galaxies very close to the 
limit C max ~ Cthrcsh- 

(2) Complete redshift surveys (CETH; Colless et al. 1993a) show that about 45%±10% 
of all field galaxies in the magnitude range of our target sample (21.25 < B < 22) 
have redshifts between 0.16 and 0.37. Our detection rate, 24 out of 49, is consistent 
with this. Given the Poisson error in the number of galaxies (20%) and large 
scale clustering of faint galaxies (e.g., the redshift "spikes" seen in the survey by 
Broadhurst et al. 1988), we conclude that we have detected at least 75% of the 
target galaxies within our surveyed redshift range, and possibly all of them. 

Hence we believe that we are measuring LWs for a well-defined, representative sample 
of distant field galaxies. 

2.5 Measuring the [Oil] Linewidths 

In measuring linewidths for the galaxies in our sample, we must account for several 
factors: 

• We do not know the shape of the line profile a priori. This is because the shape 
of the rotation curve is unknown, because the emission from ionized gas is lumpy 
and often asymmetric with respect to the dynamical center of the host galaxy, and 
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because not all of it comes from the flat part of the rotation curve (cf. Schommer 
et al. 1993; also Sec. 3.1). 

• The [On] (3727A) line is a doublet whose line flux ratio i?[oii] depends on the 
electron density in the line emitting gas. This ratio must be determined from the 
data. 

• The rest frame separation of the [On] doublet [(AA) = 2.75 A] is not much larger 
than the instrumental resolution and/or the typical kinematic line broadening 
observed in our galaxy sample. 

We cope with these complications by choosing a specific model for the line profile of 
the doublet: 

[^-^] 2 i , D „ r [A-A ze ] 2 
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with X z = (l + z)\ , X Z e = (l + ^)(l + e)A , a z = {l + z)<J v , a ze = (1 + z)(l + e)a v , and 
Ao = 3726.05 A. This form incorporates the known doublet separation expressed as a 
fractional difference, e = (AA)o/Ao = 7.3805 x 10~ 4 , and ensures identical rest frame 
velocity dispersions, a v (in A units), for the two doublet components. It approximates 
the instrumental line profile by a Gaussian whose width, = 0.82 A, has been 

determined from the widths of Cu-Ar and Fe-Ar arc lamp lines. 

By choosing the functional form given in Eq. (1), we define the galaxy's LW as the 
a of the best fitting Gaussian to the observed emission line profile. The parameters to 
be fit are the line amplitude Jo, the [On] doublet flux ratio R\om, the redshift z, and 
the rest frame intrinsic dispersion of the emission line a v (a v [A] = a v \km s _1 ]A /c). 
Further, a linear continuum baseline is fitted simultaneously over the 80 pixels neigh- 
bouring the emission line. 

The best fit parameters (Jq, -R[oii]> z , °\>) are determined by minimizing x 2 , 



2 



X 2 = f ~^ mod ^ - Ik ] (2) 

pixel k V A./fc J 

directly in pixel space (cf. Rix & White 1992). The quantities I k and AI k are the 
observed intensity and noise levels, respectively, at the k-th pixel with wavelength A. 
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The "la" error bars on the best fit values of the individual parameters are determined 
by finding the region of parameter space for which x 2 ~ Xmm = 1- The upper and 
lower error bars are calculated separately. The two error estimates are comparable for 
most parameters and for most objects, except for low values of a v (see upper panel 
of Fig. 2) and we tabulate the geometric mean of the two error estimates, A° bs . 

The best fit line parameters for all 24 galaxies are given in Table 1. The corre- 
sponding fits are superimposed onto the data in Figure 1; Figure 2 shows two of the 
fits in detail. 

The line luminosities, L[oii], and equivalent widths, EW, listed in Table 1 are 
not based on flux calibrations. Instead, they have been calculated from the exposure 
time, assuming an overall instrumental efficiency of 0.9%. Non-photometric conditions 
during the observations and inaccuracies in the fibre throughput calibration result 
in a factor of two uncertainty in the determination of L[oii] and EW. Note that 
for positive definite quantities such as the dispersion, the "best fit" value is not an 
unbiased estimator of the true value when the error is comparable to the measured 
value (cf. Wardle & Kronenberg 1974). In this regime, the probability distribution 
can be approximated by a mean of o~ v = y/o* t fit - (A° bs ) 2 and a dispersion of A° bs , 
where A° bs is the geometric mean of the upper and lower errors derived from the x 2 
fit. The best fit values of the velocity dispersion, o~ Vj at, are presented in Table 1 and 
in Figure 4; In the analysis of Sec. 3, however, we use the corrected values. Note that 
errorbars, A° bs , listed in Table 1, Figure 4 and Equation 5 are the geometric mean 
of the upper and lower errorbars. 

The main observational results based on these fits are displayed in Figure 4, 
which shows several properties of the 24 galaxies with detectable [On] emission as a 
function of their LW. This plot shows that the emission lines of most galaxies are 
resolved, but that all LWs are small: a v < 100 km s -1 . There are weak correlations 
between the observed LW (a v in Fig. 4) and other global galaxy properties. Firstly, 
the LWs appear to be somewhat larger for galaxies of greater [On] emission line 
luminosity. Secondly, the median LW tends to increase with increasing continuum 
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(B band) luminosity. The statistical significance of the trends and a comparison with 
corresponding measurements in the local gal.axy population will be given in Sec. 3 
and 4. 

The [On] doublet flux ratio, -R[oii], for most galaxies is found to be near the low 
density value of 1.45 (for n e <> 10 25 cm" 3 ; see Osterbrock 1989, p. 134), consistent 
with the line ratios found in Hn regions. Five of the 24 galaxies in our sample have 
doublet line flux ratios, R[ou] « 1.32, indicating gas densities significantly (> 2a) 
in excess of 100 cm -3 . These objects are shown by the open symbols in Figure 4. 
The [On] flux in these may arise predominantly from an active nucleus, as the LWs 
and inferred gas densities are comparable to those of LINERS (Ho, Filippenko and 
Sargent, 1993). In this small, but non-negligible fraction of our blue galaxy sample 
(5/24), the large emission line EW may not indicate vigorous star formation, but 
rather an AGN. This is in agreement with the spectroscopic results of Tresse et al. 
(1996) who find that the integrated emission line flux ratios for ~ 20% of blue field 
galaxies are inconsistent with standard Hn region spectra, and are similar to the 
ratios found in LINERs. Since we are interested in what the [On] linewidths might 
tell us about the global gas kinematics on kiloparsec scales, we exclude these five 
objects from the subsequent analysis. 

3 ANALYSIS 

In this section, we describe a set of steps that are essential for understanding what our 
LW measurements might be telling us about luminosity evolution in galaxies. First, 
we show that the line profiles must be broadened by global orbital motion (rotation) 
of the ionized gas in the galaxy; this is done by showing that the observed LWs are 
much larger than expected from the turbulent motions in individual star forming 
complexes. Second, we quantify the bias in the measurement of the rotation speed v c 
from the LW measurement by using the observed properties of nearby small galaxies, 
the potential local counterparts of the distant galaxies in our sample. For a null- 
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hypothesis prediction we use Fabry-Perot (FP) datacubes of these nearby galaxies to 
simulate the expected distribution of measured LWs, a v at given v c for these galaxies 
if they were observed at z ~ 0.25 with our AAT/AUTOFIB fibre setup. Third, we 
devise a maximum likelihood test to determine what range of linewidth-luminosity- 
colour relations are consistent with the data. 

3.1 Do the Measured Linewidths Reflect Turbulent Gas Motions? 

3.1.1 Relation between Linewidth and Emission Line Luminosity 

Melnick et al. (1989) have shown that there is a well determined relation between the 
FL3 line luminosity, Lnp, and the mean (Gaussian) LW, Wl, for giant Hn regions and 
"Hn galaxies". Presumably, these LWs arise from turbulent motions within the Hn 
regions, powered by star formation activity and the resulting supernovae. Melnick et 
al. find that this relation takes the form: 

log(oZ) = log(oSl) + rj \og{L n p/ L 40 . 5 ) (3) 

where 040.5 is the mean LW at L-^p = L40.5 = 10 4a5 ergs s _1 cm~ 2 and the slope of 
the relation is rj = 0.2. This relation holds over the range 38.5 < log(L H/ g) < 42, with 
a remarkably low scatter of only A L ~ 0.06 in log(o- L ). The emission line broadening 
due to turbulent motion is expected to be very similar for different emission lines. 
To check whether our L\om-a v data for distant galaxies are consistent with the local 
Lnp-o~ v relationship, we use the observed mean line flux ratio for field galaxies at 
(z) ~ 0.25 (Tresse et al. 1996): (L[oii]/I/h/3) ~ 4.5, with a scatter of a factor of two. 
Applying this mean line flux ratio and adjusting the Melnick et al. (1989) fit to H = 
70km s _1 Mpc _1 , we expect log(cr 40 .5) = 1.28 for the [On] line — i.e., a L ~ 19km s _1 
for a galaxy with an [On] emission line luminosity of L[oii] = 10 40 ' 5 ergs s~ x cm" 2 . 
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3.1.2 Likelihood Analysis 

We apply a maximum likelihood test to our galaxy sample in order to determine the 
best fit parameters [log^o.s), rj, A L ] and their associated errors. The relation given 
in Eq. (3), with an intrinsic scatter in log((7i), predicts a LW distribution: 

where cfj, is the mean velocity dispersion predicted for an "Hn galaxy" with [On] line 
luminosity L[oii]- On the other hand, the probability distribution of the (true) LW, 
a v , given a measured value of <r° bs and its error, A° bs , is approximated by 

v^(a ) - 1 exp ( ^~< hS y ) ( 5 ) 

For any parameter set [log(a 40 .5), rj, A L ], the probability of making the observation 
[log(L [0 ii]), < bs ], is given by 

Ppog^iom), < hs \ = I v ohs {<y)vT d W)da . (6) 

The likelihood of a L[ n]-o\, relation, given N observations, is then 



N 



log (5351), 77, A L =£logp(log(L [0 „]), < bs ) . (7) 



i=i 



The best parameters are determined by maximizing the likelihood in Eq. (13) and 
their confidence intervals are derived from the distribution of 2(£ max — C) (cf. Wilks 
1962), which is asymptotically distributed as a x 2 distribution of n degrees of freedom 
and where n is the number of parameters considered simultaneously. The maximum 
likelihood analysis yields the following results: (a) a well determined zero point of 
the relation, log(o"4 .5) — 1-62 ± 0.05 or <r 40 5 = 42 ± 3 km s -1 , which is larger by 
7a than the value in the Melnick et al. (1989) relation [log(<7 4a5 ) = 1.28 or = 
19 km s -1 ]; (b) a slope rj = 0.19 ± 0.08, (c) a scatter A L = O.ll&io, indicating that 
most of the observed scatter is attributable to measurement error. Figure 5 shows 
the data points, the best fit (bold solid line) and its uncertainties (dot dashed lines), 
and the Melnick et al. relation (dotted line). 
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A direct likelihood ratio test shows that our best fit relation is inconsistent with 
the relation found by Melnick et al. (1989) at the > 99.99% level. At an [On] line 
luminosity of L[oii] = 3 x 10 40 ergs s" 1 cm" 2 , the mean LW observed in our distant 
galaxy sample is a factor of two higher than that predicted by the <7„-L H /3 relation 
for local Hn galaxies. Put differently, local Hn galaxies, whose turbulent LWs are 
comparable to those observed for the distant galaxies should have 50 times higher 
line luminosities. This suggests that the observed [On] LWs in our distant galaxy 
sample are not determined by the turbulent motions within star forming regions, 
but rather reflect the global rotational motion of a disk of ionized gas, as they do in 
samples of local spiral galaxies. 

3.2 Relating the Measured Linewidth to the Galaxy Rotation Speed 

We calculate what distribution of measured linewidths, p(o~ v \v c ), we should expect 
for typical local galaxies of a given rotation velocity, v c , if we observed them at 
z = 0.25 with AAT/AUTOFIB. Such a mapping must account for the lack of spatial 
information in fibre spectroscopy and the effects of: 

• A finite fibre aperture (D — 2 "A) and fibre pointing errors of ~ '.'3 (a conservative 
assumption as engineering tests indicate a 0'.'2 pointing error). 

• The shape of the rotation curve [v c (R) ^ constant] 

• Clumpy and asymmetric spatial distribution of the ionized gas emissivity 

• The absence of information about the disk inclination, % 

• Seeing (FWHM ~ 1 ".!) 

• Fitting a Gaussian model line profile even though the "true" [On] line profiles of 
galaxies may have a variety of different shapes 

3.2.1 Construction of Galaxy Kinematical Models 

As a first step, we construct a complete, two-dimensional model for the gas kinemat- 
ics of each of three well-resolved local galaxies, ESO 215-G39, ESO 437-G34, and 
ESO 323-G73. The model specifies the mean velocity at every radius R and the line 
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emissivity at every point (R, 6) in the disk plane. "Reduced" FP data for the three 
local galaxies were kindly provided by T. Williams and P. Palunas. The basic prop- 
erties of these galaxies are listed in Table 2. The galaxies were chosen because they 
cover a range in absolute magnitude and rest-frame colour similar to the members 
of the distant sample. The reduced FP data for each galaxy consist of: (a) the Ha 
"flux image", I(x, y), as a function of position in the plane of the sky, and (b) the 
velocity field, v pTO j(x, y), at points of sufficient line flux [I(x,y) ~ at many points 
in the image]. The velocity field data are too sparse to allow the usual tilted ring 
fit; we instead resort to a more restricted model, one in which the galaxy is assumed 
to be axisymmetric and coplanar, and to have a rotation curve which can be well 
approximated by the parametric form: 



where x = R/R . Despite the awkward appearance of the above functional form, each 
parameter has a simple interpretation: v c is the velocity scale, Ro is the "turnover" 
radius (or core radius), 7 determines the sharpness of the turnover (the higher the 
value of 7, the sharper the turnover), and the index (5 specifies the power-law behavior 
of the curve at large radii. This form of the rotation curve has a linear rise at small 
radii with a slope: dv rot /dR ~ v c / Rq. In practice, the rotation curves are nearly flat 
at large radii (\(3\ < 0.1), and we set f3 — in our model fits in order to have a well 
defined velocity scale. 

In addition to the four rotation curve parameters, we fit the galaxy's center 
(xq, yo), mean recession velocity, inclination, and major axis position angle. All 
nine parameters describing the velocity field are well constrained by the data for 
each of the three local galaxies, and these can be used to deproject the Ha image, 
I(x, y) — > I(R, 9). Some of these parameters are listed in Table 2. The best fit rota- 
tion curve and azimuthally averaged Ha flux distribution are shown in Figure 6. This 
figure also shows the size of a 2'/l fibre aperture if these local galaxies were observed 




(8) 
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at z — 0.25 (D = 7.6 kpc). The aperture extends to twice the turnover radius Rq for 
two of the galaxies, while -Rfibrc ~ Ro for the third. 

3.2.2 Projection of Kinematical Models 

The deprojected models of the three local galaxies, consisting of the ionized flux im- 
age, I(R, 9) and the rotation curve, v TOt (R), are projected to simulate the appearance 
of a z = 0.25 galaxy viewed from any arbitrary direction. To calculate the probability 
distribution, p(a v \v c ), of measuring a width a v in a galaxy with a characteristic rota- 
tion speed v c , we need to specify the distribution of disk inclination angles, p[cos (i)]. 
A magnitude limited sample of dust free disk galaxies is expected to have a uniform 
distribution in cos(i). 

For nearby spiral galaxies, however, it is well established that an inclination- 
dependent "internal extinction correction" must be made in calculating the true 
luminosity from the observed brightness, to account for dust within the disk of the 
galaxy. The magnitude of this correction is still under some debate (cf. Rix 1995), 
but the B band correction given in the Third Reference Catalog (de Vaucouleurs et 
al. 1991, hereafter RC3) is likely to be an upper limit: Ab{i) = 1.5 min[log 10 (a/6), 1], 
where the axis ratio b/a equals cos(i) for an infinitely thin circular disk. Galaxies at a 
given inclination % with "observed" (extincted) luminosities in the interval [L, L + dL] 
have "true" (unextincted) luminosities in the interval [e As ^L, e As ^(L + dL)]. Since 
the value of A B (i) increases for more edge-on disks, a magnitude limited sample may 
be skewed towards face-on galaxies. For sub-L^ galaxies, however, this effect is small 
if the faint-end power law exponent of the luminosity function, a, is close to —1. If 
the luminosity function is a power- law with a = — 1 there are equally many galaxies 
in the interval [L, L + dL] and [yL, '-/(L + dL)}. Specifically, if we adopt the faint-end 
slope for the local galaxy luminosity function a ~ —1.2 (Efstathiou et al. 1988), the 
inclination distribution is only slightly skewed, with (cos(i) = 0.54), compared to 
(cos(i) = 0.5) for the dust-free case. Further, the average velocity reduction due to 
projection, (sin(i)) differs by only 6% from the dust-free case. 
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The presence of dust also implies that a given apparent magnitude corresponds to 
a higher intrinsic luminosity than if the galaxy were free of dust. The average value 
of the internal extinction is (Ab) 0.4 for the skewed cos(?) distribution (compared 
to (A B ) = 0.6 for a uniform distribution), and we adopt this value in the subsequent 
comparison of our distant galaxy sample to local samples (Sec. 4). 

We simulate the AAT/AUTOFIB fibre observations by drawing an inclination 
from p[cos (i)] and a random major axis position angle <y? ma j and by calculating the 
projected velocity f pr oj, at each point in the sky plane (x, y) of the projected model of 
a local galaxy. To account for turbulent dispersion intrinsic to individual Hn regions, 
the line profile at each point is assumed to be a Gaussian with a velocity dispersion 
of 15 km s _1 . 

Each pixel is given weight in proportion to its Ha emission line flux and the 
seeing is simulated by convolving each flux point (spatially) with a circular Gaussian 
of FWHM = 1 ".7. The emissivity weighted line profile is then obtained by integrating 
over the area of the (projected) galaxy model corresponding to the 2'.'1 (diameter) 
aperture of the AUTOFIB fibres. This aperture corresponds to a radius of 3.8 kpc 
at a redshift of 0.25. The effect of pointing error is simulated by shifting the center 
of the aperture (over which the integration is carried out) in a random direction, 
with the size of the shift drawn from a Gaussian distribution with rms dispersion 
0'/3. The integrated line profile for a given pair of viewing angles (i, <p ma j) is fitted 
to a Gaussian in a least-squares sense; this is a non-trivial step if the true profile 
shapes are not Gaussian. Summing the projected galaxy model (of given v c ) over 
all orientation angles (i, y? ma j) produces the probability, Pmc^^c), of measuring a 
dispersion cr v with our fibre setup. 

The probability distribution Pmc{&v\vc) differs for the three galaxies because of 
differences in the shape of the rotation curve, in the spatial distribution of the emission 
line flux from ionized gas, and, most importantly, in the rotation velocity scale v c . 
Also, we need to predict Pmc{.&v\v c ) for any v c in a certain range, even though we only 
have a small number of calibrating galaxies. We do this by assuming that the shape 
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of Pmc is constant and that the distribution scales with v c . Mathematically speaking, 
we assume that Pmc(o~ v \v c ) = p(cr v /v c ). The probability distributions for the three 
local galaxies are combined and the resulting mean p(cr/v c ) is shown in Figure 7 along 
with the rms spread between the three curves. Note that this wide range of observed 
LWs, a v , at a given v c , would arise even for infinite signal-to-noise observations with 
our technique. 

The shape of p{o~/v c ) deserves some comment, as it has a broad distribution with 
a mean value of o~ v /v c ~ 0.6. These features arise from a combination of factors. 
For an edge-on galaxy with a uniform flux distribution in the plane of the disk, a 
perfectly centered aperture that is large enough to include the flat part of the rotation 
curve yields an integrated emission line profile is roughly rectangular; the Gaussian 
that best fits a rectangle extending between — v c and +v c has a width a v ~ 0.75 v c . 
This explains why p(<J v /v c ) falls off for cr v /v c J> 0.8. Any asymmetry in the shape 
of the integrated line profile, due to either a non-uniform flux distribution or fibre 
centering error, reduces further the fitted value of a v . The finite fibre size and central 
concentration of the ionized flux distribution tend to weight the rising part of the 
rotation curve at small radii at the expense of flat part at larger radii, and this too 
reduces the measured a v relative to v c . The effect of seeing tends to increase the 
value of <j v /v c since it causes flux from the outer parts of the galaxy (at extreme 
velocities) to spill into the fibre aperture and flux from the central part of the galaxy 
to spill out of it. The tail at the low end of the distribution (a v /v c < 0.4) arises 
from face-on galaxies. For a randomly oriented galaxy sample, the mean reduction in 
velocity due to projection effects (compared to an edge-on sample) is (sin(i)) ~ 0.85. 
All the above factors — the definition of a v , line profile asymmetries (caused by spatial 
inhomogeneities in the ionized flux and by fibre centering errors), central weighting 
of the galaxy (due to finite fibre size and centrally peaked ionized flux distribution), 
smearing due to seeing, and inclination effects — are of comparable importance, and 
all of them except seeing will reduce a v /v c . 

While we measure [On] for the distant galaxies, the p(cr v /v c ) was derived from the 
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Ha distribution in local galaxies. It is interesting to ask how the spatial distributions 
of [On] and Ha compare in local galaxies, since that determines the applicability of 
the derived p(a v /v c ) to the distant galaxy sample. In a study of the star forming 
regions in 14 nearby spiral galaxies, Zaritsky et al. (1994) find that the [On] flux 
is less centrally concentrated than the Ha flux. This is caused by a radial gradient 
in the intrinsic Ha/ [On] fluxes of Hn regions (which, in turn, is related to their 
radial metallicity gradient) and in the amount of (selective) dust extinction (both 
increase with increasing radius). Both these effects cause our estimate of a v /v c to 
be smaller than the "true" value, making our best fit value of the rotation speed, v c 
(see Sec. 3.3.2) somewhat of an overestimate. This would imply that the discrepancy 
between the mean rotation speed of our distant galaxy sample and that of local 
galaxies is larger than we indicate in Sec. 4. 

3.3 The Linewidth Luminosity Colour Relation 

3.3.1 Fitting Procedure 

We are now in a position to test whether our observations are consistent with 
any linewidth-luminosity (LWL) relation, or more generally, with any linewidth- 
luminosity-colour (LWLC) relation of the form: 

log[^ c (M B , B-R)] = log [v c (-19, 1)] -r)(M B + 19) + ([(B-R)-1] , (9) 

where v c (— 19, 1) is the mean circular velocity of a galaxy with absolute blue mag- 
nitude Mb = —19 and colour B — R = 1, corresponding to our sample median. 
The dependence of the mean circular velocity on galaxy luminosity and colour are 
specified by the slopes, r\ and (, respectively. 

If the LWLC relation has an intrinsic (logarithmic) scatter of in velocity, then 
the probability of finding v c , for a given luminosity and colour, is 

dv c 1 / (log[v c ]-log[v c (M B ,B-R)])\ dVc 
p[v c \M B , B-R] — = exp — * — — '— — .(10) 



c 
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As long as A^ <C 1, it matters little whether the scatter is specified in the velocity 
or in its logarithm. The predicted probability of measuring a LW o~ v for a galaxy of 
absolute magnitude M B and colour B — R is then 

roc 

p p * cd (a v \M B , B - R) = / mc((Tv\vc) p[v c \M B , B - R] — ^ , (11) 

Jo v c 

where the averaging and bias introduced by our observing technique (see Sec. 3.2) 
are incorporated through Pmc(°~v\ v c)- As in Sec. 3.1.2, this distribution must then 
be convolved with the measurement uncertainties in order to yield the probability of 
each galaxy observation, a° bs , given an assumed LWLC: 

roo 

pK bs |M B , B - R) = / p oh \<j v ) f™ A {<j v \M B , B-R) da , (12) 

Jo 

where p obs (<7„) is defined in Eq. (5). The likelihood, £, of such a LWLC given all N 
observations is then 

N 

£M-19, 1), 77, C, A„] = ]T \og[p(a: bs \M B , B-R)] . (13) 

i=i 

3.3.2 Likelihood Limits 

The most likely values for the above parameters are determined by maximizing the 
likelihood in Eq. (13); their confidence intervals are derived from the distribution of 
2(£ m ax — £) (cf. Wilks 1962). Figure 8 summarizes the result by showing two projec- 
tions of the likelihood contours, that include the best fit value, in the [v c (— 19, 1), 77] 
and [v c (— 19, 1), Q planes, illustrating the dependence of v c on M B and (B — R), 
respectively. The analysis yields best fit values of: v c (— 19, 1) = 66 ± 8 km s _1 , 
7] = 0.07 ± 0.08, and C = 0.28 ± 0.25. The projections of the best fit onto the 
M B — i] and M B — ( plane are compared to the data in Figure 8. The likelihood limits 
of the parameters in the M B — rj and M B — ( planes are shown in Figure 9. 

The "zero point" of the LWLC relation is well determined. Note that the internal 
dust extinction is included in these models and enters in two ways: (1) it leads to 
a 0.4 mag correction for the mean 5-band internal extinction for inclined spirals 
(Sec. 3.2.2), and (2) it results in a slight skewing of the inclination distribution towards 
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face-on galaxies, changing the inferred linewidth by ~6%; these two effects "cancel" 
each other partially. There is some covariance between the zero point and the slope, 
i] = — d\og(v c ) /dMs, because the median absolute magnitude of the distant galaxy 
sample is somewhat brighter than our reference value, M B = — 19. 

As the top panel of Figure 9 shows, a wide range of slopes, rj = — d\og(v c )/dM B , 
is consistent with the data, including rj — 0, the case of no luminosity dependence of 
v c . This is caused in part by the limited absolute magnitude range of our sample of 
galaxies (AM# ~ 2 mag) and in part by the breadth of the p(cr v /v c ) distribution. Note 
that the typical LW-luminosity slope observed for local samples, rj ~ 0.12 (Bothun 
et al. 1985; Fukugita et at. 1991), is also consistent with our data. 

The bottom panel of Figure 9 shows the colour dependence of v c through the 
likelihood contours in the [i> c (— 19, 1), (} plane, at fixed rj. The dependence of v c 
on B — R, suggested by the bottom panel of Figure 4, is not significant: ( = 
d\og(v c )/d(B -R)= 0.28 ± 0.25 (68% confidence limits). 

In order to test the sensitivity of these results to the exact shape of the proba- 
bility distribution, p(a v /v c ), we have repeated the above likelihood analysis with the 
p(cr v /v c ) distributions derived separately for each of the three local calibrator galax- 
ies. The three probability distributions are roughly similar in shape as shown by the 
error bars in Figure 7. The resulting 68% likelihood regions for any given parameter 
overlap substantially for these three realizations, indicating that the results for the 
three are consistent with one another. 

The intrinsic scatter in the LWLC relation of distant galaxies [defined in Eq. (10)] 
is only weakly constrained: A„ < 0.32 (90% upper limit). This may be understood as 
follows. The probability distribution, p(cr v /v c ), for any given v c is very broad even if 
A„ = 0. In the presence of finite intrinsic scatter in the LWLC relation, the resulting 
distribution of a v is a convolution of the intrinsic v c spread (Gaussian with dispersion 
A„) with p(a v /v c ). The results are largely independent of the intrinsic scatter A„, 
as long as it is at least a factor of two smaller than the spread in observed LWs, 
p(cr v /v c ), introduced by our LW measuring technique. This implies that the results 
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(including the size of the confidence regions) are insensitive to the value of A^; we 
fix this parameter at its best fit value, A„ = 0.15, in the rest of the analysis. Note, 
the intrinsic scatter is also consistent with zero. 

4 COMPARING THE LINEWIDTH-LUMINOSITY RELATION FOR 
DISTANT VERSUS LOCAL GALAXIES 

We now turn to a quantitative comparison of the LWLC in our distant sample to the 
properties of present day galaxies. This requires LW measurements in a set of local 
galaxies whose B luminosities and colours are similar to those of the distant field 
galaxies in our sample. Unfortunately, most galaxy samples targeted in Tully-Fisher 
studies are designed for optimal distance estimates, but do not provide an unbiased 
estimate of the LWLC relation for a statistically well defined population. 

4.1 Relating the HI Linewidth W 2 o to the Optical Linewidth v c 

A large number of local galaxy LW measurements are based on Hi single-dish data, 
usually defined to be the width of the Hi line profile at 20% of the peak intensity, 
W20 (or equivalent measures such as W50). To relate W20 to the optically measured 
"circular velocity" v c , it is customary to set: 

v c = 0.5fW 20 /sm(i) . (14) 

There are two main reasons why the correction factor / differs from unity. First, 
turbulent motions broaden the Hi profile. Second, rotation curves are not perfectly flat 
at large radii, and the Ha-emitting ionized gas and neutral Hi gas (used to measure v c 
and W 2 o, respectively) may sample different radii. The relative importance of these 
two effects is still under debate and we resort to published Ha and Hi studies to 
obtain an empirical calibration of the v c - W20 relation. One of the largest samples of 
optical rotation curves is available from the survey by Mathewson et al. (1992a). In 
their analysis Mathewson et al. (1992b) they find that v c (Ha) ~ 0.94 (0.5 W 50 /sin(i)) 
for galaxies with M B ~ —19; Rubin et a/.(1985) and Courteau (1992) find that on 
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average (v c ) = v( W50) within a few percent, albeit for more luminous galaxies. If 
these results are combined with W50 ~ 0.9 W20 (e. (/.Bothun et a/.1985) this implies a 
correction factor of / = 0.86 ± 0.03. Schommer et al. (1993) applied the turbulence 
correction advocated by Tully & Fouque (1985) to a sample of galaxies for which they 
had both optical (Ha) and Hi measurements (corresponding to / = 0.83), but found 
that a smaller correction (/ = 0.9) was needed to bring the two measurements into 
agreement. In summary, these studies indicate that a correction of / = 0.86 ± 0.04 
should be applied for galaxies with Mb ~ — 19 to go form H / 2 o/sin(i) to 2v c . 

In the subsequent comparison with local data, we adopt a correction factor of 
/ = 0.86, which corresponds to correcting W 2 o downward by 14%. 

4.2 Existing Tully-Fisher Studies 

We consider two B band Tully-Fisher samples of nearby disk galaxies — those of 
Bothun et al. (1985) and Fukugita et al. (1991) — for comparison to our distant galaxy 
sample. The member galaxies of both of these samples are distant enough to yield 
sufficiently good relative distances relative to our (z) ~ 0.25 sample, once scaled 
to the same Hubble constant. Further, both these Tully-Fisher samples contain a 
sizeable number of galaxies with luminosities as low as M B ~ —18, and thus overlap 
in luminosity (and colour) with our sample of distant galaxies. 

The B band study of Coma cluster galaxies by Fukugita et al. (1991) yields a mean 
rotation velocity, v c (— 19) = 102 km s" 1 , for galaxies with M B = —19. Bothun et al. 
(1985) list B — V colours for their galaxies, allowing us to mimic the bj—r F < 1.2 color 
cut for the distant field galaxy sample. The resulting mean colour of this subsample 
is B — V = 0.36, comparable, after K-corrections and color transformations (Fukugita 
et al. 1995) to the mean (observed) colour of our (z) ~ 0.25 galaxy sample. The zero 
point of the LWL relation for the colour-restricted Bothun et al. galaxy subset is not 
significantly different from that of the whole sample, v c (— 19, 1) = 102 km s" 1 , and 
is identical to Fukugita et aVs 1991 value. 

As the the LWs in these local Tully-Fisher samples were derived from Hi measure- 
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ments, W20, the quoted value for the "optical" rotation speeds, v c (— 19) = 102 km s -1 , 
include the 14% downward correction from Section 4.1. The parameters of the best 
fit LWL relations for these two local Tully-Fisher samples are shown in Figure 9, 
with a correction for W 20 — > v c . 

4.3 Local Line-width Luminosity Colour Relation from RC3 

Most nearby Tully-Fisher samples, including Fukugita et al. (1991) and Bothun et al. 
(1985), are restricted to a narrow range of Hubble types and usually exclude irregular 
galaxies. Our (z) ~ 0.25 sample, however, has no morphological restriction. Since the 
mean v c at a given M B may vary with Hubble type (Rubin et al. 1985), the local 
Tully-Fisher galaxies may not be directly comparable to the distant ones. 

As a remedy, we have selected from the RC3 catalog a set of local galaxies with 
measured W 2 o, ignoring their Hubble type, requiring only: 

a. Axis ratios b/a < 0.7, so that the uncertainty in the LW sin _1 (i) deprojection 
factor is small 

b. Recession velocities cz > 1000 km s" 1 and angular separations of > 10° from the 
center of the Virgo cluster, so that our distance estimate, D = cz/Hq, is not 
severely affected by peculiar velocities 

The resulting RC3 sample includes galaxies as faint as M B = — 17 and many galaxies 
with — 18 > M B > —20, the absolute magnitude range spanned by our (z) ~ 0.25 
sample. 

We construct a LWL relation for the blue (B — V < 0.5, (B - V) = 0.41) 
and red (B — V > 0.5, (B — V) = 0.66) subsample separately. This colour cut 
corresponds to the colour threshold, bj — r F < 1.2 used to select the distant galaxy 
sample. Both blue and red RC3 samples show a clear LWL trend, despite a scatter 
of almost a factor of 2 in v c at fixed M B . This scatter appears to be dominated by 
disk inclination errors. The zero point of the best fit LWL relation for the blue RC3 
sample, v c (— 19, blue) = 110 km s^ 1 , is in reasonable agreement with the zero points 
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of the morphologically-restricted Fukugita et al. (1991) and Bothun et al. (1985) 
samples. As discussed above, these zero point estimates include a 14% downward 
correction from 0.5 W^ 20 /sin(?) — > v c . The best fit slope of the RC3 blue galaxy LWL 
relation, r\ = 0.11, is also consistent with that of the Fukugita et al. and Bothun et 
al. samples and with that of the distant galaxy sample. 

The mean rotation velocity of M B = —19 galaxies in the red RC3 sample, 
v c (— 19, red) = 129 km s _1 , is significantly higher than in the blue sample, although 
the LWL slopes of the two samples are identical. This trend in the local v c versus 
colour relation corresponds to a slope, ( = +0.11 (after converting the rest-frame 
B — V colour for local galaxies to a redshifted bj — r F colour at z — 0.25). 

The colour dependence of v c in local galaxies has an interesting implication. If 
the difference between the mean rotation speeds, v c (— 19, 1), of distant and local 
galaxies is caused by galaxies being brighter in the past, one also expects them to 
be somewhat bluer in the past (due to a younger mix of stars). In other words, the 
appropriate local counterpart of a z = 0.25 galaxy with Mb = —19 and B — R = 1, is 
not only one that is fainter but also one that is redder (B — R> 1). Since the mean v c 
increases for redder galaxies locally, the offset in v c between distant and local galaxies 
(and the amount of luminosity evolution it implies) would be even larger than we 
calculated in Sec. 4.2 had we taken colour evolution into account. 

In summary, comparison of our linewidth data on blue field galaxies at (z) ~ 0.25 
with nearby galaxy data leads us to conclude that, at the > 99% confidence level, 
the distant galaxy sample does not have rotation speeds as high as that of local 
samples with similar photometric properties (B luminosity, B — R colour). Taken at 
face value, the best fit rotation speed of blue M B = — 19 galaxies, v c (— 19, 1), is 35% 
smaller (at least 25% smaller at 90% confidence) than the local value. It is likely 
that this difference can be attributed to an epoch of increased star-formation, which 
lead to a larger luminosity at a given LW. Assuming the local B-band slope for the 
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LWL relation (rj pa 0.12), the off-set in v c {— 19, 1) correponds to a magnitude offset 
of AMb = 1.5 mag, with a lower limit of 0.8 mag (99%). 

5 CONCLUSIONS AND FUTURE WORK 

We have presented results from an exploratory project to determine the relation 
between the kinematic and photometric properties of galaxies at cosmologically sig- 
nificant distances and to compare it to the relation found for local samples. The target 
of our study is a statistical sample of blue, sub-L* field galaxies in the redshift range 

0. 16 < z < 0.37, for which we have obtained high dispersion spectra to determine 
their [On] emission LWs, integrated over an aperture of ~ 7.6 kpc diameter. The [On] 
LWs are resolved (er ^ 30 km s _1 ), yet the measured dispersions a v are < 100 km s _1 
for all sample galaxies. For several reasons we believe that we have detected the 
[On] emission and have measured LWs for most galaxies within the magnitude, color 
and redshift range; the sample is nearly complete. With the imposed colour cuts 
the sample constitutes the blue half of the galaxy population in this luminosity and 
redshift range. 

The primary goal of our project is to test the null hypothesis that these distant 
galaxies have the same characteristic circular velocities, v c , as local galaxies of the 
same absolute B magnitude and colour. We have used Ha Fabry-Perot data on local 
galaxies to mimic our fibre observations, and to thereby calibrate the biases/errors 
associated with the measurement of v c from emission LWs. The following main con- 
clusions arise from this analysis: 

1. At a given emission line luminosity, the observed [On] linewidths of distant galax- 
ies are too large to be compatible with the relation between line luminosity and 
(turbulent) LW observed for nearby giant Hn regions and Hn galaxies. Further, 
80% of the our sample members have [On] doublet flux ratios that indicate gas 
densities of <J 100 cm -3 , consistent with Hn region spectra. Hence, the LWs of 
most objects reflect the orbital motion of ionized gas in the galaxy's potential 



28 Rix, Guhathakurta, Colless, Ing 

well. In the remaining 20% of objects, the emission line flux may arise from an 
active nucleus. 

2. We have tested whether distant blue galaxies within a narrow range of absolute 
blue luminosity have the same mean v c as local galaxies of the same absolute 
magnitude and colour. In order to make a realistic comparison, we have simulated 
what emission LWs would be measured at z ~ 0.25 for local galaxies of a given v c , 
assuming that the shape of the rotation curve and spatial distribution of ionized 
gas are comparable in distant and local galaxies. We find that, at a given lumi- 
nosity, Mb ~ —19, the distant galaxies have rotation speeds that are about 35% 
smaller than expected from photometrically identical, local samples. The LWs of 
the two samples are inconsistent at the > 99.9% confidence level. 

3. The most likely explanation for the LW offset is that galaxies have undergone 
considerable luminosity evolution. Assuming the local B-band slope for the LWL 
relation (r) 0.12), the off-set in v c (— 19, 1) correponds to a magnitude offset 
of AM B = 1.5 mag, with a lower limit of 0.8 mag (99%). If spectral evolution 
accompanies this luminosity evolution, the magnitude offset may be even larger. 
However, at this point it cannot yet be ruled out that the LW differences in the 
samples arise, at least in part, from a very different v c — >[On] mapping. Such 
differences could arise from either the shape of their rotation curve, and/or the 
spatial distribution of the line emitting gas. 

This systematic investigation of low-luminosity, blue field galaxies, complements 
corresponding studies of massive, red cluster galaxies. Franx (1993) and van Dokkum 
and Franx (1996) studied the fundamental plane of early type galaxies in two clusters 
at z — 0.19 and z = 0.4. They found that galaxies of a given mass scale were brighter 
in the past by an amount that is consistent with a formation at high z and subsequent 
"passive evolution". Also, Lilly et a/.(1996) find that the luminosity function evolves 
most strongly for low-luminosity blue galaxies. The present data indicate that this 
evolution of the luminosity function can be traced to the mass-to-light ratio evolution 
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of individual galaxies. However, we infer values of v c for galaxies of Mb = —19 which 
are not as low as predicted for some dwarf starburst scenarios (e.g. Babul and Rees, 
1992; Babul and Ferguson, 1996). 

The present study represents only one step towards understanding the internal kine- 
matics of distant field galaxies. Follow-up should address the following questions: 

a. Is there evern more direct evidence that emission LWs reflect ordered rotation of 
ionized gas within these galaxies? We are currently pursuing this question with the 
help of Fabry-Perot imaging datacubes for a small sample of distant field galaxies 
(Ing et al. 1996; see also Vogt et al. 1993, 1996). 

h. What is the spatial distribution of the emission line flux in these distant galaxies? 
Is it the same as for local galaxies, or are the observed small LWs in part due to 
a much more centrally concentrated flux distribution? A direct answer to these 
questions would require high angular resolution, narrow band images tuned to the 
redshifted [On] line. 

c. Can the scatter in the linewidth-luminosity relation be reduced if we correct for 
galaxy inclinations? A reduction in scatter is expected if, as we argue, the LW 
is dominated by gas motions in a disk. Inclinations of sufficient accuracy can be 
obtained from HST images. 

d. What is the slope of the linewidth-luminosity relation for distant galaxies when 
we include a wider range of absolute magnitudes? We are analyzing high disper- 
sion spectroscopic data obtained with the AAT/AUTOFIB instrument of nearly 
100 galaxies spanning a broader range of Mb and colour than the sample presented 
here. 

e. Is there a similar offset in the LWL relation when the luminosity is measured in a 
different bandpass? The amount of luminosity evolution is expected to be smaller 
in the rest-frame I band than in B. Near infrared (if -band) photometry is being 
carried out to test this hypothesis. 
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Table 1. Summary of the Observational Results 
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The error listed is the geometric mean of the upper and lower error bar. 

The EW has been estimated from the observed [Oil] line flux, the instrument efficiency and the bj magnitude of 



the galaxy. 

c The five galaxies where the line ratio is significantly (2<r) below 1.32 have been labelled with an asterisk. Their line 
ratios 

indicate gas densities in excess of 100cm -3 , and may arise from an AGN, rather than HI regions in the galaxies' 



disks. 
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Table 2. Local Galaxies for Establishing p(a\v c ) 



Name i>[km/s] M B B — R v circ i? [kpc] i [°] 7 

ESO 215-G39 4335 -19.8 0.8 153 1.88 47 2.2 

ESO 323-G73 5016 -19.9 1.0 157 1.83 48 2.3 

ESO 437-G43 3801 -17.7 0.8 104 3.50 53 3.7 



